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Are the spin axes of stars randomly aligned within a 
cluster? 
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ABSTRACT 

We investigate to what extent the spin axes of stars in young open clusters are aUgned. 
Assuming that the spin vectors he uniformly within a conical section, with an opening 
half-angle between A = 0° (perfectly aligned) and A = 90° (completely random), we 
describe a Monte-Carlo modelling technique that returns a probability density for this 
opening angle given a set of measured sinz values, where i is the unknown inclination 
angle between a stellar spin vector and the line of sight. Using simulations we demon- 
strate that although azimuthal information is lost, it is easily possible to discriminate 
between strongly aligned spin axes and a random distribution, providing that the 
mean spin-axis inclination lies outside the range 45°~75°. We apply the technique to 
G- and K-type stars in the young Pleiades and Alpha Per clusters. The sini values are 
derived using rotation periods and projected equatorial velocities, combined with radii 
estimated from the cluster distances and a surface brightness/colour relationship. For 
both clusters we find no evidence for spin- axis alignment: A = 90° is the most prob- 
able model and A > 40° with 90 per cent confidence. Assuming a random spin-axis 
alignment, we re-determine the distances to both clusters, obtaining 133 ± 7pc for 
the Pleiades and 182 ± 11 pc for Alpha Per. If the assumption of random spin-axis 
alignment is discarded however, whilst the distance estimate remains unchanged, it 
has an additional percent uncertainty. 
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1 INTRODUCTION 

Most authors considering the statistics of orbital or rota- 
tional stellar motion assume, where relevant, that angular 
momentum vectors are randomly orientated. It is possible 
however that the physical processes of star formation lead 
to a preferred axis of rotation over the scale of an individual 
star forming region (SFR). This might arise if the direction 
of average angular momentum of the molecular cloud giving 
rise to the SFR has a significant influence on the resulting 
angular momentum of individual stars - for instance, if gas 
were constrained to collapse along strong, large-scale mag- 
netic fields threading the cloud. 

Historically, little has been discussed either theoreti- 
cally or observationally about the possibility of spin-axis 
alignment during star formation. This would require a rel- 
atively undisturbed collapse along magnetic field lines with 
little disruption from turbulence or dynamical interactions 
(e.g. Shu, Adams & Lizano 1987). Using circumstellar disc 
orientation as a proxy, some studies have suggested pref- 
erential spin alignment with the ambient magnetic field in 



SFRs (e.g. Tamura & Sato 1989; Vink et al. 2005), but oth- 
ers have found no evidence for disc axis alignment (Menard 
& Duchene 2004). 

A second important reason for assessing the degree of 
spin-axis alignment is that measurements of projected equa- 
torial rotation velocities (wsini, where i is the unknown in- 
clination of the spin-axis to the line of sight) and rotation 
periods can be combined to provide a powerful method to de- 
termine the distances (Hendry, O'Dell and Collier-Cameron 
1993; Jeffries 2007a; Baxter et al. 2009), radii (Jackson, Jef- 
fries & Maxted 2009) or star formation histories (Jefi'ries 
2007b) of young clusters. Such statistical analyses must as- 
sume that the orientation of spin axes are random. If the spin 
axes in an individual cluster were in fact partially aligned, 
this would change the intrinsic sin i distribution for the clus- 
ter, producing biased estimates of distances, radii and age 
spreads. 

In this paper we use published rotation data for the 
young Pleiades and Alpha Per open clusters to investigate 
to what extent spin axes may be aligned once the star for- 
mation process has finished. In Section 2 we discuss how well 
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spin-axis orientation can be determined using measured ro- 
tation periods and projected equatorial velocities. In Section 
3 we present a parameterised model for the observed sini 
distribution resulting from a group of stars with partially 
aligned spin axes and show how well such a model can be 
used to determine the degree of alignment from simulated 
datasets. In Section 4 we compare the models with measured 
sin i distributions for G- and K-type stars in the Pleiades and 
Alpha Per clusters. Section 5 presents and discusses the re- 
sults of our analyses, including new, independent estimates 
of the distances to these clusters under the assumption that 
the spin axes are randomly aligned, and in Section 6 we give 
our conclusions. 



2 OBSERVATION OF SPIN AXES 

Using current measurement techniques it is not possible 
to observe stellar spin-axis orientation directly. For young, 
magnetically active, spotted late-type stars it is possible to 
measure their period of rotation, P, from rotational mod- 
ulation of their light curves and their projected equatorial 
velocity, usin i, from spectral line broadening. These can be 
used together with photometric data and an independent 
measure of distance to determine sini, the sine of the angle 
between the observer and the spin axis 
P 

sini = — (usini) , (1) 

2-kR ^ ^ ' ^ ' 

where the stellar radius R can be estimated from the surface 
brightness and distance. Such observations give no informa- 
tion on the azimuthal direction of the spin axis, only an 
estimate of the inclination, i, which is also degenerate be- 
tween i and n — i. 

The effect of restricting observations to measurements 
of sini is illustrated in Fig. 1. The left-hand panels indi- 
cate the intrinsic distribution of spin-axis vectors. The cen- 
tral panels show what would be observed given the lack 
of azimuthal information and the degeneracy in i discussed 
above. The right hand panels show the cumulative sini dis- 
tributions that would be observed (see section 3.1). 

In the first case, where there is no preferred orientation, 
the uniform distribution of spin axes would be observed as 
a uniform inclination distribution over a hemisphere from 
to 7r/2. The next case is more complicated. Here, the spin 
axes of a group of stars are distributed uniformly over a con- 
ical region about a central cone axis with defined inclination 
and azimuthal direction. This would be observed as a larger 
circular region on the surface of the hemisphere that is sym- 
metric about the line of sight; the azimuthal information 
being lost. 

The third and fourth cases show how conical regions 
representing the same degree of alignment would be ob- 
served at increasing average inclinations. The lack of az- 
imuthal information causes spin axes represented by a rel- 
atively small cone area to sweep out a large area of the 
measured hemisphere. This concept of swept areas is ap- 
proximate. A detailed treatment is given below, where the 
variation in number density over the swept area is modelled 
to calculate the cumulative probability density as a function 
of the opening angle of the cone. A, and its mean inclination, 
a. The point here is to emphasise that azimuthal and inclina- 
tion angle degeneracies conspire, especially when convolved 



with measurement errors (see next section), to hamper the 
recovery of the underlying sini distribution. 



3 MODELLING THE sin/ DISTRIBUTION 

We will use the term "sin i distribution" for a group of stars 
to express the set of numbers specifying the angle between 
the rotation axis and the observer's line of sight. Two such 
distributions can be considered: 

(i) a true sini distribution which depends only on the 
distribution of the spin axes of a set of stars, 

(ii) a measured sini distribution which depends both on 
the true distribution and the uncertainties and limits that 
apply to measurements of the base parameters used to de- 
termine sini. 

3.1 The true sini distribution 

The true sini distribution depends on the distribution of 
the stellar spin axes over the celestial sphere. The simplest 
case is a uniform (random) distribution in which case the 
cumulative probability distribution depends on the area of 
the celestial sphere between an angle and i 

-Ptruc = 1 ~ cos i for i = to 7r/2 , (2) 

where the superscript o denotes the case of a uniform dis- 
tribution. 

A simple way to represent an aligned distribution is to 
assume that spin vectors are uniformly distributed over a 
conical solid angle and zero elsewhere (see Fig. 1). The cone 
angle. A, which corresponds to half the opening angle, de- 
termines the degree of alignment. A small cone angle means 
stars have nearly parallel spin. A large cone angle (A ~ 7r/2) 
corresponds to a uniform distribution. The mean inclination 
of the stars within the cone is represented by a. 

The equivalent sin i distribution can then be calculated 
using a Monte Carlo method as follows. A set of spin axes 
are specified at angles On and <;/>„, relative to the cone axis. 
The cumulative probability distribution depends on the an- 
gle between the cone axis and spin axis, 9„, as 

P(e') = (1 - cos6i)/(l - cosA) for 61 = to A (3) 

and random values of 9„ are generated as 

0„ = cos"'(l-7?„(l-cosA)) (4) 

where R„ is a random number between and 1. The prob- 
ability distribution of the angle around the cone axis is 
uniform, allowing random values of (f) to be generated as 
0,1 = 2ttR'„ (where _R„ and R'„ are different random num- 
bers). The inclination relative to the line of sight is cal- 
culated by considering the triangle formed by unit vectors 
along the spin axis and the line of sight with respect to the 
cone axis. As the line of sight is at an angle a with respect 
to the cone axis, then 

cos in = sin a sin On cos + cos a cos On ■ (5) 

Since measurements of inclination derived from projected 
rotational velocities cannot distinguish between i and n — i, 
the effective value of sin i is given by 

sinin = sin(cos~"'^(j sin a sin 6„ cos (p„ + cos a cos 0„\)) . (6) 
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(a) Actual distribution of spin axes (b) Observed distribution (c) Cumulative probability 




sin i 

Figure 1. A schematic diagram siiowing (a) the distribution of spin vector orientations, (b) the corresponding distribution that would 
be deduced from sin i measurements and (c) the resultant cumulative probability distribution of measured sin i (solid line) compared to 
that found for a uniform distribution of spin axes (dashed line) . 



A Monte Carlo method is used to determine the ex- 
pected distribution of sini, for a given A and a. A set of 
sini values is evaluated for random values of 6 and (f>. The 
results are then ordered to define the cumulative distribu- 
tion function of sini. This cumulative distribution can then 
be used to determine a representative set of sin i values by 
generating random numbers between and 1. 

The right hand panels in Fig. 1 show the cumulative 
sini probability distributions calculated for a uniform dis- 
tribution and for well aligned distributions with A — 25° and 
various values of a. For A = 25° the true sini distribution 
is significantly different from that of a uniform distribution. 
Thus with sufficient sin i measurements it should be possible 
to differentiate between well-aligned and uniform distribu- 
tions, irrespective of the inclination a. 

Things become less clear when A is increased to, say, 
60°. In this case the results at low (a ~ 20°) and high 
(a ~ 75°) inclinations could still be differentiated from a 
uniform distribution. However, at intermediate inclinations 
the sinz distribution becomes quite similar to that for a 
uniform distribution. Thus a weakly aligned spin-axis dis- 
tribution with large A might only be discernible if its mean 
inclination is either low or high. 

Of course, these simple considerations have so far ig- 
nored the alterations to the observed sin i distribution that 
are imposed by selection effects in the data and by measure- 
ment uncertainties, which are discussed in the next section. 



3.2 Measured sin i distribution 

An expression for the observed value of sin i is obtained from 
equation 1 and by assuming that the stellar radius is propor- 
tional to the product of its angular diameter and distance. 

sin lobs = fc-Pobs("sini)obs/(^obs-Dcst) , (7) 
where fe is a constant (appropriate to the units used) of 

^ = ^truc-f-^truc/(-ftruct'truc) , (8) 

Pobs is the observed period of a star, (i;sini)obs is its mea- 
sured projected equatorial velocity, Aohs its angular diame- 
ter (derived from a magnitude and colour via a Barnes-Evans 
relation - see section and Dost is the distance to the star 
estimated by some independent method. 

Uncertainties in the observed vsini, period, angular di- 
ameter and distance estimate can be represented as Gaus- 
sian distributions with normalised standard deviations of Sp, 
Sv, 5a and Sd, such that Pobs ~ -Ptruc(l + SpU) etc. where 
[/ is a random numbers drawn from a Gaussian distribution 
with mean of zero and unit standard deviation. Hence we 
can write 

fl + 5pvUi\ 

\l + OADU2 ) 

where 8pv = Sj, + 6^, Sad = \/ S\ + Sj, 

and Ui and U2 are different random numbers with a mean 

of zero and a standard deviation of 1. For the moment Spv 
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Figure 2. Simulations of the cumulative probability distributions of observed sini assuming normalised uncertainties Spy = Sad = 0.1 
and a distribution of spin-axis orientations that is either a uniform (random) or is aligned within a cone angle A = 30°. The rows show 
results for mean inclinations of a = 30°, 60° or 80°. The columns show the effects of altering sininjin- The dotted line shows the effect of 
setting (Dsini)niin = 0. Note the value of sini can exceed unity due to uncertainties in the measured parameters used to determine sini. 



and Sad can be considered as empirically derived constants. 
Their values are discussed in section 

In addition to measurement uncertainties there are 
thresholds below which either period or wsini cannot be 
measured. These thresholds are equivalent to a lower limit 
to sin i below which rotational modulation and periodicity 
would not be detected, and a resolution limit defining a 
threshold for v sin i detection. The latter is reasonably well 
defined from the observational data, but there are significant 
astrophysical uncertainties in the former - e.g. the latitude 
distribution of spots on a young, active star (see Jeffries 
2007a for a discussion). We chose to represent this observa- 
tional bias as a simple cut-off value sinimin, below which a 
period could not be obtained for a star. In Figure 2 we show 
that the value of sinimin has a non-negligible effect on the 
observed sini distribution, so it is treated as a free param- 
eter in the following analysis and allowed to vary between 
zero and 0.71. That is to say we do not specify a lower limit 
of sini, and at worst we expect to be able to measure the 
period and v sin i of stars with inclinations 45° and above. In 
section 5.1 it is shown that this is justified by the available 
observations. 

To model the effects of the resolution limit for projected 
radial velocity measurement, (ii sin i)niin, we require an esti- 
mate of_th£distributi^ of utruc. The approach used here fol- 
lows Ijeffriesi |20073) whereby the intrinsic Utrue distribution 
is represented as a combination of a uniform distribution 



and an exponential decay. For the Monte Carlo analysis a 
fraction 7 of velocities are drawn from a uniform distribution 
between zero and Umax and the remainder from a cumula- 
tive exponential distribution of the form P(v) = exp(— 
In practice the sini distribution is not sensitive to the ex- 
act form of the Vtrue distribution so parameters defining the 
velocity distribution {'^, P^Vmax and vsmimin) can be esti- 
mated by matching the distribution of (iisini)obs. 

To take into account the threshold values of v sin i and 
sini in our simulations, any realisation with sin itrue < 
sinimin or (v sin i)obs < (wsini)min is excluded from the sini 
distribution, since these would not be present in an observed 
data set. 

Figure 2 shows Monte Carlo simulations of the cumu- 
lative sin i distribution with typical levels of uncertainty: 10 
per cent in combined period and projected radial velocity 
and 10 per cent in combined angular diameter and distance 
(see section 4.6). Results are shown for three values of mean 
inclination, a — 30°, 60° or 80° and for sinimin = 0.2, 0.4 or 
0.6. The solid lines show results for A = 30° and the dashed 
line shows the probability density for a uniform distribution 
of spin axes (A = 90°). Introducing uncertainties reduces the 
difference in sin i distribution between the A = 30° cone and 
the uniform distribution (compare with Figure 1). However, 
the distributions are still quite different for either low or 
high values of a. For intermediate values (45° ^ a ^ 75°) 
it may still be possible to resolve the difference provided 
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(a) Surface plot 




(b) Contour plot 
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Figure 3. Analysis of the sini distribution measured on a group of 36 stars to determine the probability density of parameters describing 
the underlying distribution of spin axes as a function of cone spread (or solid angle) (1 — cos(A)) , mean inclination, (1 — cos(«)) and 
threshold in sin j. The measured distribution was derived for stars in Pleiades (see Table 3 and section 4.4) 



sinimin ~ 0.2. Above this, the effect of increasing sinimin is 
to compress the distributions along the sin i axis effectively 
erasing the distinction between the aligned and random spin- 
axis distributions. 

These results were calculated for a velocity distribution 
described by 7 = 0.33, /3 — 37kms~^, Vmax = 140kms~^ 
and (i;sini)jnin = 3.6 kms~^ (appropriate for the Pleiades - 
see section 4.5). Also shown as a dotted line in Figure 2 are 
results calculated assuming no lower cut off in vsmi. This 
produces only small changes in the modelled sini distribu- 
tions, justifying the use of a simple representation for the 
velocity distribution. 



3.3 Fitting parameters to measured distributions 

fn our models there are three unknown parameters that de- 
fine the observed sin i distribution - the opening half-angle 
of the cone. A, that describes the degree of spin-axis align- 
ment, the mean inclination a and sin imin which we will refer 
to as r. The next step in the analysis is to determine how 
the sini distribution of a group of stars can be analysed to 
determine the underlying parameters (A,Q,r), of which we 
are most interested in determining A. 

Elements of the measured sin i distribution are ordered 
to produce a cumulative distribution function. This is then 
compared with a Monte Carlo model sin i distribution using 
a Kolmogorov-Smirnov (K-S) test. This gives an estimate of 
the probability that the measured data could be drawn from 
the same distribution as the model data set. The larger the 
estimated probability the more likely it is that the model 
parameters (A, a, r) represent the measured data. 

The K-S probability is calculated for all values of 
Ai, a-,, Tfc and normalised by setting the integral over all pa- 
rameter space to be unity, to give the probability matrix 
A'"(Ai, Oj, Tfc). To visualise results it is useful to sum over 
one or two of the independent variables to produce contour 
plots of probability density for two parameters of interest or 
line plots for one parameter of interest respectively. Figure 3 
shows typical results for a set of 36 stars in the Pleiades (see 
section 4). Figure 3a shows a surface plot of the probability 
density as a function of "cone spread" , (or cone solid angle - 



defined as 1 — cos A), and the sini threshold, r. These results 
show that there is a range of A and r that give reasonable 
fits to the measured data. Very roughly, they indicate that 
the cone spread lies between «0.4 and 1 (corresponding to 
55° < A ^ 90°) and that r < 0.5. The modelling suggests 
then that the measured data are most likely drawn from a 
nearly uniform distribution, but that there is still a finite 
probability of a partially aligned distribution of spin axes. 

To put this on a quantitative basis we can look at a 
a contour plot of cone spread against 1 — cos a (Fig. 3b). 
The contours contain the labelled percentage of the summed, 
normalised K-S probability. This shows that if the cone in- 
clination is small or close to unity (a close to 0° or 90°) 
then the results could only be consistent with a large cone 
spread 1 — cos A > 0.6, or A > 66°). However, for interme- 
diate values of inclination the measurements are consistent 
with almost any value of cone spread. Figure 2 shows why 
this is. At intermediate inclinations there is a much smaller 
difference between the sin i distributions for well-aligned and 
randomly orientated spin axes. 

Figure 3c shows a summary of the results for the mea- 
sured sin i distribution. The solid line shows the probability 
density of cone spread (1 — cos A), fntegrating under this 
curve gives a 90 per cent probability that the cone solid 
angle is greater than 0.21, corresponding to a cone angle 
A > 38°. We can also say with 90 per cent confidence that 
the threshold in sini is less than 0.5. 



3.4 Conditions for discrimination of cone angle 

The example in Figure 3 demonstrates how an observed sin i 
distribution could be analysed to investigate the underlying 
distribution of spin-axis orientation for a group of stars in a 
cluster. This section considers how useful the method might 
be in practice. Specifically it considers: 

(i) Under what conditions will this method correctly re- 
cover the underlying distribution of spin-axis orientation? 

(ii) What effect does sample size and uncertainty in in- 
dividual sini measurements have on the accuracy of the 
method? 
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Figure 4. The probability density as a function of cone spread (1 — cos A) estimated from the analysis of the sini distributions of 300 
simulated stars. The curves show the results for simulated datasets with partially aligned rotation axes (A = 30°) with different mean 
inclination angles, compared with a simulated dataset with random spin-axis orientation (A = 90°, dashed line). The left hand plot show 
the comparison for cone inclinations, ol = 15°, 30° and 45°. The right hand plot shows the comparison for a = 55°, 65° and 75°. 



To address the first question the Monte Carlo method is 
used to generate sini distributions for two scenarios. The 
first represents the case where there is significant align- 
ment of spin axes of stars in the cluster corresponding to 
A = 30°. The second case represents a uniform distribution 
with A = 90°. The sini values are produced for sets of 300 
stars with typical levels of uncertainty in (sini)obs, namely 
5pv ~ Sad = 10 per cent, a representative utruc distribu- 
tion and T = 0.4. The cumulative sini distributions are then 
analysed using the method described above to determine the 
probability density of cone spread, which should reflect the 
underlying distribution of spin axes. 

The results are shown in Fig.4 for increasing values of 
a. For low values of a (< 45°) the distribution is peaked 
towards low values of 1 — cos A, roughly corresponding to the 
input value of A. At high values of a 75°) this is also the 
case. However at intermediate values, the probability density 
becomes more uniform indicating that the sin i distribution 
could have arisen from almost any value of A. The reasons for 
this were discussed in section 3.1 and need not be repeated 
here. 

Where the probability density falls to zero at one end 
or other of the distribution, we can calculate limits on A. For 
example, integration under the dashed line in Fig. 4 (where 
the input value of A was 90°) returns a 90 per cent confl- 
dence limit that A > 50° . This limit is the flgure of merit we 
choose to characterise a sin i distribution, indicating whether 
it results from an aligned distribution or not. 

To determine the effect of sample size, measurement 
uncertainties and the threshold r, a range of Monte Carlo 
simulations were performed. These were divided into runs 
with A — 90° and A = 30° . The results of analysing the sin i 
distributions are listed in Tables 1 and 2. 

The main conclusions from these simulations are: 

(i) For an input A — 90° (Table 1), the lower limit that 
can be placed on the recovered A value is quite insensitive to 



Table 1. Lower limit to the cone angle, A, found from input 
simulated datasets with random spin-axis orientation. 



No. stars 


Threshold in 


90 per cent lower limit 


analysed 


sini, T 


of cone anj 


jle, A (deg) 


Uncertainties, 


Spv & Sad 


0.10 & 0.10 


0.15 & 0.15 


10 


sini > 0.2 


> 30±6 


> 29±7 




sini > 0.4 


> 30±4 


> 29±5 


30 


sini > 0.2 


> 39±4 


> 33±9 




sini > 0.4 


> 36±6 


> 33±5 


100 


sini > 0.2 


> 46±6 


> 40±7 




sini > 0.4 


> 47±2 


> 48±3 


300 


sini > 0.2 


> 47±3 


> 48±5 




sini > 0.4 


> 49±5 


> 44±5 



the sample size. Observing many hundreds of stars does not 
give much improvement over sample sizes of ~ 30. Neither 
are the results sensitive to the actual value of r or the exact 
value of the measurement uncertainties. 

(ii) For a strongly aligned distribution (input A — 30°, 
Table 2) we find that upper limits to A are only obtained if 
the mean inclination a ^ 45° or a ^ 75°. Again, the gains 
to be made by observing very large samples of stars are not 
very significant, although samples of ~ 100 may be required 
to identify a strongly aligned spin-axis distribution when a 
is large. 



4 MEASURED sin J DISTRIBUTIONS 

In this section we use the previously developed ideas to 
model the observed sini distributions obtained from pub- 
lished measurements for stars in the young Pleiades and Al- 
pha Per open clusters. These clusters have ages of ~ 120 Myr 
and ~ 80 Myr respectively, so many of the G- and K-stars 
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are rotating fast enough to produce measurable rotational 
broadening of their spectral lines - giving (t)sini)obs - and 
produce rotational modulation caused by magnetic starspots 

- giving Pobs- Angular diameters are determined from re- 
ported V and K magnitudes using a surface brightness re- 
lation calibrated by interferometry (Kervella et al. 2004). 

4.1 Measured data for the Pleiades and 
Alpha-Per clusters 

Even in the well-studied Pleiades and Alpha-Per clusters, 
there are surprisingly few stars where v sin i and P are both 
known. The database for galactic open clusters (Mermilliod 
1995) currently lists 296 stars in the Pleiades cluster with 
known v sin i, and 56 with known P, not all of which overlap. 
Similarly, there are 247 stars in Alpha-Per with known v sin i 
and 66 with known P. Table 3 shows the (usini)obs and 
Pobs for 44 stars in the Pleiades and 38 in Alpha Per that 
have been identified as cluster members and for which this 
simultaneous information is available. 

Table 3 also lists the available stellar photometry, which 
is subsequently used to check for unresolved binarity and es- 
timate the angular diameters. We tabulate a mean V mag- 
nitude and the (peak-to-peak) amplitude of the light curve 
modulation used to find the stellar rotation period. The ap- 
parent K magnitudes are taken from the 2MASS catalogue 
(Cutri et al. 2003). A small ofi'set is applied to convert these 
to the CIT photometric system used in the evaluation of an- 
gular diarneter ( see Section |4]2} Kcit ~ K2AIASS + 0.024, 
ICarpenteil (|200ll ). 

4.2 Estimation of angular diameter 

The method used to estimate stellar angular diameter in 
equation 9 follows that of O'Dell, Hendry & Collier Cameron 
(1994). A Barnes-Evans relationship (e.g. Barnes & Evans 
1976) determines the angular diameter of a star from its 
measured apparent magnitude and a colour index. Whilst 
the approach of O'Dell et al. is appropriate for the current 
application, the calibration data they used for the Barnes- 
Evans relationship has been superseded. In addition, O'Dell 
et al. used the B — V colour index, but this is now known to 
be a systematically unreliable temperature indicator in mag- 
netically active, spotted cool stars (e.g. Stauffer et al. 2003). 
We prefer to use V — K, which is potentially more precise 
and appears less affected by stellar activity and metallicity 

- the latter being a potential source of uncertainty in any 
calibration sample. 

Kervella et al. (2004) provide a recalibration of the 
Barnes-Evans relationship based on angular diameters of 
main sequence and sub-giant stars measured by interferom- 
etry. They give the following relationship between angular 
diameter. A, (in units of arcseconds) de-reddened magni- 
tude, Ko and colour {V — K)o 

\og{A) = (0.5170 ±0.0017) + (0.0755 ± 0.0008)^1^ - 7^)0 
-O.2K0 . (10) 

This is valid for dwarfs of spectral type A0-M2 and has an 



Table 2. Upper limit to the cone angle, A found from input sim- 



ulated datasots 


with cone angle A 


= 30° and 


r = 0.4 


No. stars 


Mean 


90 per cent 


upper limit 


analysed 


Inclination 


to cone angle, A (deg) 


Uncertainties, 


Spv & Sad 


0.1 & 0.1 


0.15 & 0.15 


10 


15° 


< 41±4 


< 52±8 




30° 


< 64±12 


< 65±8 




45° 


< 74±6 


< 76±10 




75° 


not resolved 


30 


15° 


< 39±6 


< 42±4 




30° 


< 53±6 


< 54±3 




45° 


< 72±6 


< 73±11 




75° 


< 74±6 


< 80±3 


100 


15° 


< 32±1 


< 34±5 




30° 


< 51±4 


< 58±4 




45° 


< 64±3 


< 68±2 




75° 


< 66±14 


< 56±16 


300 


15° 


< 30±3 


< 30±5 




30° 


< 49±1 


< 47±6 




45° 


< 66±2 


< 67±3 




75° 


< 53±17 


< 53±11 



exceptionally small intrinsic dispersion (< 1%), although in 
our case there are other sources of uncertainty in the es- 
timated angular diameter, notably uncertainty in apparent 
magnitude and interstellar extinction. To allow for interstel- 
lar extinction the corrections given by Rieke and Lebofsky 
(1985) (Av = 3.1E{B - V) and E{V - K) = 2.7AE{B - V)) 
are applied to the above expression to give 

log(yl) = 0.5170 + 0.0755(1/ - K) - 0.2K 

-0.136£(B- V), (11) 

where V is the peak apparent magnitude for the variable 
star {Vmean — AV/2) (see Table 3). The rationale here is 
to use the brightest value of V to represent an unspotted 
photosphere. A mean value of colour excess, E(B-V) is used 
for each cluster. A value of 0.032 is taken from An, Terndrup 
& Pinsonne ault(2007) for the Pleiade s and a value of 0.10 
taken from iPinsonneau lt et all ^199^ ) for Alpha Per. The 
expression for log(A) is quite insensitive to the magnitude 
of the colour excess, a very conservative uncertainty of 0.05 
in Eb-v would only change an estimated angular diameter 
by 1.5 per cent. 

4.3 Effects of binarity 

The stars in Tables 3 will inevitably include some unresolved 
binaries. As the effect of binarity will be to decrease the ap- 
parent magnitude (by up to 0.75 mag) and possibly redden 
the star with respect to the intrinsic colour and magnitude of 
a single star, equation 11 shows that unrecognised binarity 
could change the estimated angular diameter and deduced 
value of sin i. 

To mitigate this, probable binaries were eliminated by 
rejecting stars that lie significantly above the mean cluster 
sequence in a colour magnitude diagram. To do this in a 
systematic manner, reference stars identified as being clus- 
ter members in the WEBDA data base which showed 80% 
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Table 3. Estimated values of sini for stars in the Pleiades and Alpha Per clusters using; the Barnes-Evans relation in equation 11, a 
mean distance of 131.8it2.4pc and a colour excess of 0.032 for the Pleiades and a distance of 176.2±4.9pc and a colour excess of 0.10 for 
the Alpha Per cluster. The full table is available at Blackwell Synergy as supplementary material to the on-line version of this table. 



Star 


Period 


V sin i 


Apparent 


Apparent 


Colour 


Variation 


Estimated 


Estimate 


Name 


with ref. 


with ref. 


Magnitude 


Magnitude 


excess 


magnitude 


diameter 


of sin i 




(days) 


(km/s) 


Vmean 


K2MASS 


V - KciT 


AV 


A (arcsec) 




Hll 191 


3.100 a 


9.1 f 


14.38 


10.61 


3.75 


0.04 


0.047 


0.84 


Hll 253 


1.721 b 


38.2 f 


10.66 


8.95 


1.69 


0.12 


0.072 


1.27 


Hii 263 


4.820 a 


7.8 f 


11.63 


9.39 


2.22 


0.16 


0.065 


0.81 


Hii 293 


4.200 c 


5.7 f 


10.79 


9.06 


1.71 


0.02 


0.067 


0.50 


Hii 314 


1.479 d 


41.9 f 


10.56 


8.90 


1.64 


0.09 


0.073 


1.19 


Hii 320 


4.600 c 


10.8 f 


11.04 


8.87 


2.15 


0.06 


0.080 


0.87 


Hii 345 


0.723 e 


18.9 f 


11.40 


9.27 


2.11 


0.07 


0.066 


0.29 


Hii 357 


3.400 c 


10.0 f 


13.32 


10.02 


3.28 


0.07 


0.057 


0.83 


Hii 739 


0.917 e 


14.4 f 


9.44 


7.94 


1.48 


0.03 


0.108 


0.17 ** 


Hii 883 


7.200 a 


3.8 f 


13.05 


10.25 


2.78 


0.10 


0.048 


0.80 



** Stars identified as probable binaries and excluded from the analysis of the sini distribution 



The letter following the period value denotes the reference for the period and photometric data. 

The letter following the V sini value indicates its source. The sources are as follows: a j<rishimmurtlH_et_aL (_1998), 

b MariUi et al. (1997) . c Prosser et al. (19951 . d Prosser et al. (1993) . e Messina (2001) . f Queloz et al. (1998) , 

h Stauffer fc Hartmann (1987) . i Soderblom et al. (1993) , i Stauffer et al. (1989) . k Stauffer et al. (1985) 

1 Prosser et al. (1993) . m Prosser fc Grankin (1993) . n Bouvier (1996) . p O'Dell et al. (1994) 




Figure 5. (a) Selection of probable binaries from stars of known period and tisini in (a) the Pleiades and (b) Alpha Per. The solid line 
is located 0.4 magnitudes above the dashed regression line through all cluster members. Stars lying above the solid line are considered 
as probable binaries and excluded from the measured sin i distribution. 



or greater probability of membership from proper motion 
measurements were plotted in a, V versus V — K diagram 
(see Fig. 5). Stars lying more than 0.4 magnitudes above 
the linear regression line through these reference data were 
identified as probable binaries. When this criteria is applied 
to the Pleiades data in Table 3 then 8/44 of the original 
sample were eliminated. A similar analysis for Alpha Per 
eliminated 2/38 of the stars in Table 3. 



4.4 Cluster distances 

To calculate sini from equation 7 we need cluster distances. 
We adopt 131.8 ± 2.4 pc (distance modulus 5.60 mag) for 
Pleiades and 176.2 ± 4.9 pc (distance modulus 6.23 mag) 
for Alpha-Per, base d on multi-colour main sequence fit- 
ting (see Table 1 of iPinsonneault et al.l l|l998l )). There has 
been controversy over the Pleiades distance following a 
much lower parallax-based m easurement using the Hippar- 
cos satellite of 120.2 ± 1.9 pc (|van Leeuwen fc Hansen Ruij 
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I1997I : Ivan Leeuwenll2009l ). The higher distance, consistent 
with the results of main-sequence fitting and the distance to 
binary stars, has generally been been adopted in the litera- 
ture. These include: a distance modulus of 5.60 ± 0.07 mag 
from radial veloc ity and inter feromctric measurements of At - 
las (HD 23850) jZwahlen et al. 2004); Munari ct al] (|2004l ') 
used radial velocity and photometry of HD 23642 to get a 
distance modulus of 5.60 ± 0.03 mag, althou gh re-analysis of 
the same data bv ISouthworth et al.l (|2005l ') gave a distance 
modulus of 5.72 ±0.06 mag; parallax measurements on three 
stars in Pleiades made us ing the fine guidance se nsors on the 
Hubble Space Telescope (jSoderblom et al.ll2005h gave a dis- 
tance modulus of 5.65±0.02 mag. The source of the discrep- 
ancy between the Hipparcos results and other measurement 
techniques remains unresolved. The distance to the Alpha- 
Per cluster has not been investigated to the same extent as 
the Pleiades. The Hipparcos distance of 172. 4± 2.7 pc agrees 
well with the main sequence fitting distance in this case. 



4.5 Distributions of projected velocity 

For the Monte-Carlo analysis we need an estimate of the in- 
trinsic properties of the equatorial velocity distribution (see 
section 3.2). These are estimated from the observed vsini 
data as follows. The maximum equatorial velocity, Wmax, and 
minimum projected velocity, v sin imin are taken as the max- 
imum and minimum observed vsini values. The true distri- 
bution of equatorial velocities is parameterised in terms of 7 
and l3 (see section 3.2). These coefficients are chosen to max- 
imise the probability that a model v sin i distribution and the 
observed data set are drawn from a common distribution us- 
ing a K-S test. Values of 7 = 0.33 and /3 = 37kms^^ for 
Pleiades give a K-S probability of 0.8. Values of 7 = 0.16 and 
/3 = 54kms~^ for Alpha Per give a probability of 0.7. We 
re-emphasise that the exact choice of these parameters has 
very little effect on our Monte-Carlo models of the observed 
sini distribution (see section 3.2 and Fig. 2). 

4.6 Measurement uncertainties 

To model the measured sini distribution we assume that 
statistical variations in the parameters used to derive sini 
are Gaussian distributions (see equation [9]) with normalised 
uncertainties that are derived from the measured data and 
published uncertainties. The source papers give uncertain- 
ties for about half the periods in Tables 3. The average of 
this subset of measurements gives 5p = 0.03. The average 
normalised uncertainty for the velocities from Queloz el al. 
(2001) in Table 3 is 5^ = 0.10. Including uncertainties for 
data from other sources increases the average to 0.12. On 
the basis of these results the normalised uncertainty due to 
the combined effects of period and velocity. Spy, is taken to 
be either 0.10 or 0.15. 

The main source of uncertainty in estimated angular 
diameters is uncertainties in apparent magnitudes. Normally 
these should be relatively low, say ±0.05 mag for V and even 
less for K. However, there is additional uncertainty for the 
stars used here which exhibit rotational modulation. From 
Tables 3 the rms magnitude variation is AKms = 0.09 mag. 
Taking this as an upper bound for the uncertainty in both 
V and K then, from equation 16, the uncertainty in angular 



diameter is ±5 per cent. If we include the small intrinsic 
uncertainties in the Barnes-Evans relationship, < 1 per cent 
(Kervella et al. 2004) and colour excess this gives a total 
uncertainty of about 6 per cent in the estimated radii. 

The uncertainty in distance can be considered in two 
parts. First there is a random uncertainty due to the location 
of a star within the cluster. This depends on the radial depth 
of the cluster. We assume that the radial depth of the sample 
is similar to its tangential width on the sky. The Pleiades 
sample is located a mean distance of 0.8° from the cluster 
centre, corresponding to a normalised distance uncertainty 
of 0.014. For Alpha Per the sample is scatter over a larger 
area, corresponding to a normalised distance uncertainty of 
0.031 On this basis a value of Sad = 0.10 gives a conservative 
estimate of the combined effects of uncertainty in angular 
diameter and distance. 

In addition to the uncertainties applied to individual 
stars, there is an uncertainty, Snast, in the estimate of aver- 
age distance to the cluster. This produces a systematic error 
in the analysis, which becomes significant if the normalised 
uncertainty in the distance, Dcst is comparable with the un- 
certainty in the mean value of sini due to the combined 
effects of the other sources of error measured in A*' stars. 

foost ~ ViSlv + Sln)/N (12) 

Thus for A'^ = 36 stars, with Spvand Sad of 10 per cent, the 
uncertainty in estimated distance becomes significant at a 
level of 2.4 per cent. As the assumed distance uncertainties 
to the Pleiades and Alpha Per are 1.8 and 2.8 per cent re- 
spectively, then the sensitivity of our results to variations in 
this value must also be considered. 



5 RESULTS 

5.1 Probability density of cone angle for Pleiades 
and Alpha Per 

To investigate the degree of spin-axis alignment we look 
at the constraints we can place on the cone angle A. The 
Pleiades stars, filtered for binaries, were modelled using the 
method described in section 3. Normalised uncertainties of 
Spv = Sad ~ 0.10 were initially assumed. Figure 3 showed 
plots of the probability distributions corresponding to A, 
mean inclination a, and threshold sini value r. Figure 6 
now shows a similar plot of the probability density of cone 
spread ( 1 — cos A) for both Pleiades and Alpha Per, eval- 
uated for different levels of combined uncertainty in period 
and V sin i and for the cluster distances fixed at the assumed 
value plus or minus their formal error bars. 

For both clusters there is a clear increase in probability 
for larger values of cone spread (1 — cos A) and therefore 
of A itself. A qualitative comparison with the simulations 
shown in Fig. 4 demonstrate that this behaviour suggests a 
random spin-axis orientation, rather than one with a small 
A. By integrating the probability distributions we can say 
with 90 per cent confidence that A > 39° for the Pleiades 
and A > 42° for Alpha Per (see Table 5). The most likely 
model for both clusters is one with close-to random spin-axis 
alignment (A ~ 90° ) . Figure 3b shows that a lower value of A 
would have to be "disguised" in the data by an intermediate 
value of alignment, a (see the discussion in section 3.3). 
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(a) Pleiades 



2.0 ■ 




0.4 0.1" 
Cone spread (1-cos?.) 

Expected distance and 10% uncertainties 
Lower estimate of distance (10% unc.) 




r:..| O.G 0.8 1.0 

Cone spread (1-cos?i) 

Expected distance and Spy =15%, S^a=^0% 
Upper estimate of distance (10% unc.) 



Figure 6. Analysis of measured sini distributions for tfie Pleiades and Alpha Per. The solid line shows the probability density of cone 
spread (1 — cos A) that describes the underlying distribution of spin axes for normalised uncertainties of SpY = ^ad = 0.10. The dashed 
line shows the effect of increasing the normalised errors to 0.15 and 0.10. The dashed/dotted curves show the effects of varying the 
assumed distance to the cluster by plus or minus one standard deviation. The dotted line in the left hand plot shows the probability 
density of cone spread calculated assuming the a lower estimate of distance based on Hippacos data. 



Table 4. Constraints on the cone angle, A, and on the thresh- 
old value of sini, t for the Pleiades and Alpha Per for different 
assumed distances and levels of measurement uncertainty. 



Distance 


Uncertainty 


Cone Angle A 


Threshold r 


to cluster 


Spv 1 Sad 


90% limit 


90% Umit 


Pleiades 








131.8 pc 


0.10, 0.10 


> 39° 


< 0.48 


131.8 pc 


0.15, 0.10 


> 37° 


< 0.50 


129. 4pc 


0.10, 0.10 


> 43° 


< 0.48 


134.2 pc 


0.10, 0.10 


> 36° 


< 0.47 


Pleiades usin 


g Hipparcos distance 




120.2 pc 


0.10, 0.10 


> 40° 


< 0.63 


Alpha Per 








176.2 pc 


0.10, 0.10 


> 42° 


< 0.52 


176.2 pc 


0.15, 0.10 


> 37° 


< 0.57 


171.3 pc 


0.10, 0.10 


> 39° 


< 0.51 


181.1 pc 


0.10, 0.10 


> 44° 


< 0.54 



For both the Pleiades and Alpha Per we also find similar 
90 per cent upper limits to r of ~ 0.5, suggesting that it is 
not until they have inclinations below 30° that stars become 
unlikely to exhibit significant rotational modulation in their 
light curves. This supports the assumption made in section 
3.2 that T will always be less than 0.71, which corresponds 
to imin ^ 45°. Figure 6 and Table 4 also illustrate how 
robust these results are to alterations in our assumptions 
about the level of measurement uncertainty. Increasing the 
5pv to 0.15 does not change the probability distributions in 
any significant way and has only a modest impact on the 90 
per cent lower and upper limits to A and r. 



Changing the assumed distance could have more influ- 
ence. In Fig. 6 we show curves derived for assumed distances 
that are equal to the main sequence fitting distances plus or 
minus an error bar. For the Pleiades reducing the distance 
would yield a larger A lower limit, whereas increasing the 
distance hints at the possibility of some alignment though 
still consistent with a uniform distribution. We could dis- 
card the main-sequence fitting distance entirely and instead 
adopt the lower Hipparcos distance. This case is also shown 
in Fig. 6. Whilst the cone angle probability becomes a little 
more uniform there is still evidence of alignment. The re- 
sults in Table 4 show that A is still greater than 40°. We 
caution the reader that this should not be taken as evidence 
that these methods are insensitive to the assumed distance. 
As discussed in section 4.6, a distance accurate to a few per 
cent is really require for trustworthy constraints on A. 

5.2 New distance estimates to the Pleiades and 
Alpha Per 

If we adopt the assumption that the spin-axes really are 
randomly oriented then we can derive an independent dis- 
tance(0'Dell et al. 1994). This provides a first check for 
possible alignment of spin axes since if the distance derived 
assuming random orientation does not agree with the estab- 
lished distance there is a strong indication of partial align- 
ment which can be investigated using the method described 
above. Note the converse is not true since a matching dis- 
tance can also be produced by a partially aligned cluster at 
a favorable average inclination (45° to 75°). 

Combining the v sin i and period data in Table 3 gives 
an estimated distance to the Pleiades of 133 ± 7pc. This 
is in good agreement with the value found from main se- 
quence fitting of 131.8 ± 2.4 pc adopted in the last section 
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(|Pinsonneault et al.lll99d ). but siR nificantly higher th an the 
Hipparcos value of 120.2 ± 1.9 pc (jvan Leeuwenllioogl ). The 

uncertainty in our independent distance estimate is too large 
to completely rule out the Hipparcos result, but clearly fa- 
vors the more conventional main-sequence fitting measure- 
ment (assuming that spin axes truly have a random orien- 
tation!). 

For Alpha Per the data in Table 3 gives an estimated 
distance of 182 ± llpc. This agrees well with the distance 
given by main sequence fitting, 176.2 ± 4.9 pc, and within 
one standard deviation of the value derived from the new 
reduction of Hipparcos data, 172.6 ± 2.8 pc (|van LeeuwenI 
I2OO9I ). 

It is also interesting to compare these results with those 
of O'Dell et al. (1994) who used a similar technique to find 
the distance to these clusters, but with smaller samples. 
They found distances of 132 ± 10 pc and 186 ± 12 pc for 
the Pleiades and Alpha Per respectively. To some extent 
the agreement here is fortuitous. O'Dell et al. used an older 
Barnes-Evans calibration based on B — V data and did not 
exclude probable binaries. Using our data set with the sur- 
face brightness B—V relations of Kervella et al. (2004) would 
tend to increase the estimated cluster distance, whereas our 
exclusion of possible binaries reduces it. 



5.3 Biases in estimated cluster distances 

If the distribution of spin-axis orientation were not random 
then there may be a systematic error in the distances es- 
timated in the last section and it might be possible to ob- 
tain agreement with the Hipparcos distance for the Pleiades. 
In Fig. 7 we show the results of Monte-Carlo simulations 
which investigate by what factor we would overestimate the 
true distance for the cases of two aligned distributions with 
A = 15° and A — 45° respectively, if we were to analyse the 
sini distributions under the assumption of random align- 
ment. In each case we have assumed that the mean inclina- 
tion a is randomly distributed between a maximum value 
of 90° and a minimum value of 15° which corresponds to 
spin-axes which almost point towards the observer. We take 
15° as a practical minimum value because for lower values 
we simply wouldn't be able to measure any v sin i values. 

The mean x-axis values for the two distributions shown 
in Fig. 7 are both 0.99. That is, the distance we would esti- 
mate (on average) is not significantly biased by the assump- 
tion of random spin-axis alignment. However, the widths 
of these distributions imply a significant additional scat- 
ter. Clusters in which the spin-axes were highly aligned 
(A = 15°), with a mean inclination a = 15° would have 
their distances under-estimated by a factor of two. If on the 
other hand the mean inclination were to take its maximum 
value of Q = 90°, the distance would be overestimated by a 
maximum factor of about 1.2. 

If A is not 90°, then it seems reasonable to take a 68 
per cent confidence interval about the mean of the A — 15° 
distribution in Fig. 7 as an indication of the maximum 
additional systematic error in the estimated distance that 
is introduced by assuming random spin-axis orientation. 
Quantitatively, this amounts to an additional uncertainty 
of percent. 



0.3 f 



0.2 ■ 



0.1 - 



Minimum 
values 

a= 15° 



Well aligned stars - 
cone angle X = 15° 



Partially aligned stars - 
cone angle X = 45° 



Maximum 
values 




0.6 0.8 -.0 
Estimated/true distance 



Figure 7. Possible bias in the estimated distance to a cluster for 
well (A = 15°) and partially aligned (A = 45°) stars in a clus- 
ter. The histogram shows the frequency of estimates of distance 
(normalized to the true distance) assuming that the mean axis of 
alignment of the cluster is oriented randomly in space. 



6 DISCUSSION AND SUMMARY 

The results in this paper show that analysis of the measured 
sini values for groups of as few as 36 stars in a cluster can 
be used to investigate whether their spin axes are randomly 
oriented or well aligned in space. The information that can 
be gained is incomplete because (a) information is lost in the 
measurement technique, principally in the azimuthal direc- 
tion of the spin axes and (b) there are significant uncertain- 
ties in the measurements of parameters used to determine 
sini, which blur its distribution function. 

The simulations in section 3.4 and Fig. 4 indicate that 
a cluster of stars with random spin-axis orientation can be 
identified from the strong increase in probability density 
with cone angle. In some cases it is easy to distinguish this 
from a well-aligned distribution of spin axes, but it depends 
on the average inclination, a, of the aligned spin axes. If 
a ^ 45° then a well aligned distribution produces a dis- 
tinctive probability density that falls to zero for A < 90°. 
However, if 45 < a < 75° then the probability density of A 
becomes quite fiat, giving no clear indication of the under- 
lying distribution of spin axes. The simulations presented in 
section 3.4 also show that there are only modest gains to be 
made by observing very large samples of stars The exception 
would be a set of stars with a well-aligned distribution with 
a high mean inclination (a ^ 75°). In this case a relatively 
large sample (~ 100 stars) is required to clearly identify the 
high degree of alignment. 

The first application of this technique using data for 36 
stars in the Pleiades and Alpha Per is encouraging. For both 
clusters there is a clear increase in probability density with 
cone angle, entirely consistent with random spin-axis orien- 
tation. However, we cannot rule out partial alignment of the 
spin axes, but place 90 per cent lower limits of about A > 40° 
for both clusters. Values of A that are much less than 90° 
could only be disguised in the data if the mean inclination 
were between 45° and 75° . Our results also show that the 
analysis method is relatively insensitive to the exact levels of 
measurement uncertainty and small changes in the assumed 
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mean distance to the cluster. These results show that for 
the method of analysis to be effective the average distance 
to the cluster taken from independent measurements needs 
to be well defined with a normalised uncertainty of a few 
per cent or better. 

Should it turn out to be possible to assume that spin 
axes are always randomly directed within a cluster, then this 
allows measurements of average stellar distances, radii and 
age spreads. In this paper we have used this technique to 
derive new distance estimates for the Pleiades and Alpha 
Per of 133 ± 7pc and 182 ± 11 pc respectively. These values 
are in good agreement with the main sequence fitting dis- 
tances to these clusters. The Pleiades result is marginally 
higher than the Hipparcos parallax-based distance. If larger 
samples of stars could be measured, then unlike the test 
for random spin alignment, there are significant gains (im- 
proving by a factor of approximately %/iV) to be made in 
the precision of the mean sini value, which in turn deter- 
mines the statistical precision of the distance estimate. As 
there are far more v sin i measurements than known periods 
in both the Pleiades and Alpha Per, it would be a fruit- 
ful project to search for more rotation periods in both these 
clusters. However, if spin-axis orientation cannot be assumed 
random, then although the estimated distance remains un- 
biased when analysed under the assumption of randomness, 
there are significant systematic errors of up to percent, 
which render this technique ineffective as an independent 
means of estimating stellar radii or cluster distance. 

Well defined distances should soon become available for 
many more open clusters following the launch of the GAIA 
satellite. In addition periods are now being measured sys- 
tematically in nearby clusters (e.g. Aigrain et al. 2006). As 
these data become available it will be interesting to compare 
the distances derived from satellite based parallax measure- 
ment with those estimated from the measured period and 
rotational velocity. If all results agree within measurement 
errors then this would support the general assumption that 
spin-axes are randomly orientated in space over the scale 
length of a cluster. This is because it would become increas- 
ingly unlikely that the spins in all clusters were aligned, and 
had a mean inclination in the range 45°-75°. Conversely, if 
the spin axes in some clusters are well aligned there would 
be a reasonable chance (~ 55 per cent if we can assume 
that the mean inclination from cluster-to-cluster is a random 
variable) of observing a cluster with an average inclination 
of 45° or 75° in which case significant alignment could 
be identified using the method described in this paper. 
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